Heat pumping in nanomechanical systems 
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We propose using a phonon pumping mechanism to transfer heat from a cold to a hot body using 
a propagating modulation of the medium connecting the two bodies. This phonon pump can cool 
nanomechanical systems without the need for active feedback. We compute the lowest temperature 
that this refrigerator can achieve. 



Freezing out atomic motion by cooling matter to ab- 
solute zero temperature is a thought that has, for ages, 
fascinated both scientists and laymen alike. In atomic 
gases, techniques such as evaporative cooling can bring 
temperatures down to the submicrokelvin scale, allow- 
ing for the observation of quantum phenomena such as 
Bose-Einstein condensation. In solid state matter, the 
ionic motion takes the form of oscillations around equi- 
librium positions, and completely freezing the system (in 
the case of an insulator) means removing all lattice vibra- 
tions - phonons - leaving solely the quantum mechanical 
zero-point motion. 

The quest for observing quantized mechanical motion 
in macroscopic systems has incited several experimental 
groups in recent years [1]. In most cases, cooling is ob- 
tained by a feedback mechanism which involves optical 
or electronic sensors and some control system that acts 
directly on a cantilever. In this Letter, we argue that it is 
possible to cool a nanomechanical system without relying 
on feedback control. The mechanism we propose acts di- 
rectly on the acoustic phonons carrying heat in and out of 
the system without the need for monitoring its state. By 
deforming the lattice in the medium connecting the me- 
chanical system to its phonon thermal reservoir, one can 
pump heat against a temperature gradient by extracting 
out phonons. The mechanism resembles a classical cool- 
ing cycle of a thermal machine and its physical basis is 
time-reversal symmetry breaking. The pump works in 
both coherent and incoherent phonon regimes. 

Quantum coherent electron pumps have been studied 
extensively since Thouless's original proposal p|. For in- 
stance, using lateral quantum dots and quantum wires, 
charge [sl] , spin Q , and heat ^ currents can be created in 
the absence of bias by modulating adiabatically and pe- 
riodically in time two independent external parameters. 
In contrast, pumping massless bosons such as acoustic 
phonons is a much more subtle problem. For one, it 
is much harder to pump adiabatically phonons due to 
the lack of a large energy scale such as the Fermi en- 
ergy. Moreover, phonons not only obey a different wave 
equation but are also not conserved when scattered by 
external perturbations that couple linearly to the dis- 



placement field (i.e., a driving force). The result in this 
case is entropy generation in addition to pumping. 

In practice on can pump phonons with minimum heat 
generation by coupling quadratically to the displacement 
field, either by locally modulating the propagation ve- 
locity or by locally applying a pinning potential. An ex- 
treme example of a pinning perturbation, which preserves 
phonon number, is one that imposes Dirichlet boundary 
conditions to the displacement field at a given point in 
space. When such a perturbation travels along a quasi- 
onc-dimcnsional medium, it works as a linear peristaltic 
pump. Below, we show that this mechanism allows for 
cooling down the system to a minimum temperature Tmin 
which, in one-dimension, is given by the expression 
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with Qb — A^Sub/Stt'^c, where Th is the temperature of 
the hot thermal reservoir, A is the perturbation strength, 
c is the phonon velocity, and vb is the barrier speed. 

A scheme of the pumping cycle is shown in Fig. [1] 
where the nanomechanical system to be cooled is rep- 
resented by the left (cold) side. The local modulation 
in the phonon velocity or pinning potential works like a 
moving semireflective barrier to the phonons. In process 
A— >-B, the barrier is translated from the cold to the hot 
side of a cavitylike region. After it reaches the endpoint, 
another barrierlike perturbation is activated at the oppo- 
site side of the cavity (process B— >-C). Then, in C— >A', 
the first barrier is deactivated and phonons from the hot 
reservoir free expand into the cavity. The procedure is 
then repeated. 

Interesting issues arise out of this simple process of 
moving a reflective barrier (a "mirror") for phonons, in 
particular that of phonon pressure across the barrier. In- 
deed, a similar process to the one described above was 
used by Bartoli when he attempted to show the ap- 
plicability of thermodynamics to electromagnetism and 
raised the question of radiation pressure [6|, which in 
turn inspired Boltzmann in his studies of blackbody radi- 
ation [3]. The issue of phonon pressure is not trivial (and 
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more subtle than the case of photons) as phonons carry 
crystal momentum (q) but not obviously physical lin- 
ear momentum (denoted by p) . The connection between 
these two forms of momentum requires anharmonicity 
and is given by Pq — jdhq, where 7 is the Griineisen 
parameter of the lattice and d denotes the spatial dimen- 
sion This impacts the relation between pressure 
and energy density in a phonon gas; for instance, in the 
case of a single acoustic mode, the relation takes the sim- 
ple form p = ~ {dF/dV)j, = jE/V. 

We begin by discussing first the case of a fully reflective 
barrier. We can treat the problem as a gas of phonons, 
which we cycle according to Fig. [1] Notice that the bar- 
rier does not let heat pass through and the cooling is due 
to the removal of internal energy from the left-hand side, 
dumping it into the right-hand side, as explained below. 

The expansion A ^ B is adiabatic and reversible 
{ASr^l — 0, i.e., no heat exchange between left- and 
right-hand sides). Recalling the standard equation for 
massless bosons, dS = dE/T -\- E dV/VTd, we can re- 
late changes in energy to variations in volume. When 
the barrier moves to the right, the change in internal en- 
ergies on the two sides are E^ = Ef — PLVpipe/^d and 
E^ = E^ + PRVpipc/-fd, where Vpipc is the swept vol- 
ume. Then, once we insert the other barrier to get to C, 
we redraw the boundary of what L is. The volume of L 
changes by a factor {Vl — V'pipc)/VL- So in C we have 

- (1 - ^) EE, = El and E^^^ = ^ i?f , 
where the last energy is the one inside the "pipe" . Then, 
once the right barrier is removed in going C ^ A' , one 
redraws the boundary of what i? is, so E^ 
Er ^ E'^ + iJ^pc- Putting it all together, we have 
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where . . . stand for terms down by powers of Vpipo /Vl,r, 
and eR,L = Er^l/Vr^l are the intensive energy densities 
in the two sides. All the work done occurs in A B and 
is given by 
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where the leading term is insensitive to changes in the 
pressures pl,r as the volume expands. The unusual rela- 
tion between work and volume change shown in Eq. (|2c|) 
comes from the fact that, in our scheme, volume changes 
also require an increase in the number of unit cells, so 
that the lattice unit cell volume is kept constant (i.e., no 



compression) . The work required to add units cells leads 
to the 7d factor dividing the pressure difference. (No- 
tice that this factor is absent for photons, since 7 = l/d 
follows from Pq = ?iq. In the case of light, there is no 
undelying lattice system - an "ether" - that needs to be 
accounted for.) In the process A A! described above, 
all entropy increase occurs when the barrier is removed in 
going C — >■ A', and the second law of thermodynamics is 
satisfied. From this analysis, we can compute the energy 
flux out of the left reservoir per unit time of operation of 
the cycle: 
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where vb is the barrier speed. Here we use for total time 
the duration of the A ^ B stroke, assuming that the 
equilibration in the entropy production part C — > A' is 
fast compared to this time. 
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FIG. 1. Pumping cycle: A— >-B— s-C— >A' (see text for an ex- 
planation). A traveling lattice perturbation acts as a semire- 
flective barrier moving from cold to hot reservoir. The wide 
arrows indicate unimpeded heat flow. 

For our case of interest, cl = rjdT^^^/c'^, where 
r]d = 2gdlC{d + l)/[{4Trf/^T{d/2)] with C{z) and r(z) 
denoting the Riemann zeta and Gamma functions, re- 
spectively, while (7 is a degeneracy factor. Notice that 
the energy flux depends only on the intensive quantities 
for the system on the left (and thus on Tl), and not 
on any property on the right-hand side of the barrier, in 
particular its temperature. This is a straightforward con- 
sequence of the fact that the barrier is perfectly reflective, 
so one is not faced with the difficulty of flghting a thermal 
gradient between the hot and cold reservoirs. The ideal- 
ized situation, however, serves the purpose of displaying 
clearly the main principle of our cooling mechanism. 

Let us turn the discussion to the less idealized situa- 
tion when the barrier is not perfectly reflective, allow- 
ing some heat to be transmitted from the hot to the 
cold side. In this case, we intuitively expect that the 
slower we move the barrier, the more difficult it becomes 
to cool, because the energy transferred in the operation 
A ^ B ^ C ^ A' depends only on the volume swept by 
the barrier, but not on the rate (as long as the A ^ B 
stroke is done in a quasi-equilibrium situation, allowing 
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for thermal equilibration on both sides of the barrier). In 
addition, the longer we take to move the barrier to the 
right in the A B stroke, the more heat is transferred 
through the transmitting barrier (the total transfer scales 
linearly with the sweeping time). So let us now compute 
the heat flow through the moving barrier, and the con- 
ditions to attain net cooling for a semireflective barrier 
moving with speed vb- Hereafter, for simplicity, we focus 
on a purely one-dimensional case (d = 1). 

For concreteness, consider a "moving mirror" corre- 
sponding to a region in space where the atoms are cou- 
pled to an external short-range potential, which is lo- 
calized in space. The position of this pinning potential 
is modulated in time so as to make it travel at speed 
Vb, causing the reflection and transmission coefficients 
to depend on the red and blue shifted frequencies of 
the phonons coming from the two reservoirs. Acous- 
tic phonons in a one-dimensional chain, interacting with 
such a "moving mirror" potential of strength A, obey the 
following wave equation in the continuum limit: 

d^u{x, t) - c^dlu{x, t) = -\c 5{x - vb t) u{x, t), (4) 

It is simpler to work in the reference frame of the barrier, 
t' = t and x' = x — VBt, where the wave equation becomes 

{df -VBd^>f -c^d^] u{x',t') = -XcS{x')u{x',t'). 

(5) 

Let us consider plane wave solutions to Eq. ([5]) in the 
two regions, to the left of the barrier (with amplitudes 
and Bj) and to its right (with and B^): 



[x',t') = I due"^'' (A±e-'"^'/^'« + B±e*"="'/"^) , 



' (6) 

with vji = c — vb and vl = c + vb- The function u{x\t') 
and its partial time derivatives are continuous, but its 
partial space derivative is not. Integrating Eq. ([5]) be- 
tween 0~ and 0^ yields the remaining boundary con- 
dition. Matching the solutions on the two sides of the 
barrier using the boundary conditions yields 
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M±(w) 
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Using Eqs. ((7a)) and (I7b|) . the scattering matrix connect- 
ing incoming and outgoing amplitudes can be computed: 
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— i2uj/\ —1 
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Now, to determine the heat transmission and reflection 
coefficients, one needs to go back to the reference frame of 



laboratory (i.e., that of the reservoirs), where the Bose- 
Einstein occupation numbers of the phonons are known: 
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where 



{aja^^) = riLiuj), {b+h+) ^ nR{uj), and {ajb+) = 
(fe+^flj^) ~ 0, since phonons coming from different reser- 
voirs are uncorrelated. Thus, the heat current leaving 
the left reservoir is given by the expression 
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The quantity can be expressed in terms of 

the distributions nL,ji{uj) through the scattering matrix 
S{lo). After a few manipulations, we arrive at 
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The first line of Eq. ([TU]) is the thermal heat current 
/thermal from left to right in the presence of a nonmov- 
ing barrier. The second line, which we name /, 
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sults from the barrier motion and it is clearly zero when 
Vb ^ Q {vl — "vr ^ c). In the limit when the barrier 
amplitude is high, A 3> Tr^l, we obtain 
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Notice that this current is always negative if < and 
Vb > (with vr > Vl), thus, as expected, we are fighting 
this heat flux with the energy flux of Eq. . A net flux 
of energy is indeed possible if we satisfy 
which requires 
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As mentioned earlier, for a fully reflective barrier (A 
cxd), cooling can be obtained for any temperature gradi- 
ent. For a semireflective barrier, to leading order in vb /c, 
cooling requires Tl > Tmin, where Tmin is given by Eq. 
^ with Th = Tr. Notice that when Tl = Tr = T, 
the proposed mechanism also allows one to transfer heat 
between reservoirs provided that T/A < [1 / 2tt) ^J ^ / 2tt , 
independently of the barrier speed. 

A few remarks are in order. First, we note that the in- 
equality (I12p is independent on 7. In fact, anharmonicity 
is not essential for the operation of the cooling mecha- 
nism. Although anharmonicity is necessary for equilibra- 
tion to occur in a closed system, it is not so in an open 
system coupled to thermal reservoirs. For the latter, 
equilibration and thermalization takes place over time 
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scales of the order of the time required for sound waves to 
propagate back and forth through the system. Straight- 
forward numerical simulations of a harmonic linear chain 
of masses and springs coupled to a thermal reservoir at 
finite temperature show that, for practical purposes, fast 
equilibration is achieved when the barrier moves between 
reservoirs with a speed a few times smaller than c. This 
is important because anharmonic effects are very weak 
at low temperatures UJ, |l2| and should not significantly 
contribute to equilibration. Second, work is inevitably 
done when the barriers are activated and deactivated 
during the _B — > C and C ^ A' processes. However, 
during a fixed cycle, this work does not scale with the 
length of the cavity connecting the two reservoirs, while 
the amount of energy extracted from the cold reservoirs 
does. Therefore, the contribution of this work to the 
energy balance of the cooling process can be made very 
small for a sufficiently long cavity and due to this reason 
we neglected it in our estimates of the minimum cooling 
temperature T^in- Finally, although Eq. (fTOl) has been 
derived assuming coherent heat transport, Eq. ([3]) does 
not rely on quantum coherence. Hence, coherence is not 
an essential ingredient for our heat pump. 

Finally, let us discuss practical implementations. To 
produce a propagating barrier, it is better to use elec- 
tromechanical couplings rather than purely mechanical 
ones, since electronic controlled is both more precise and 
allows for faster switching times. Strongly electrostric- 
tive materials, in which changes in phonon dispersion are 
caused by an external electric field, could be used. In par- 
ticular, electrostrictive polymers such as poly-vinylidene 
fluoride (PVDF) , in which giant electrostriction has been 
observed |T§| , appear to be a promising class of materials 
for building phonon pumps. Like other one-dimensional 
systems, a single chain of PVDF has four acoustic phonon 
branches: one longitudinal, two transverse, and one twist 
mode. Being a highly ionic (or polar) polymer, PVDF 
has a permanent dipole moment per monomer unit which 
couples to the external electric field, leading to a gap in 
the acoustic twist mode dispersion. Therefore a local 
electric field can virtually block the torsion modes with 
frequencies below the gap from propagating in PVDF 
14| . which is equivalent to introducing an infinite bar- 
rier for such phonons in our scheme. 

In this particular implementation, we can understand 
more clearly other aspects of the phonon pump. For 
example, the insertion or removal of the phonon bar- 
rier corresponds to turning on or off the electric fields. 
Because the field causes a phonon gap for the torsional 
modes, if the insertion is adiabatic, the energy required 



to do so is given by the phonon energy density that is 
excluded from the barrier region. As long as the barrier 
is much narrower than the length of the channel, this 
energy can be much smaller than the energy pumped as 
the barrier is pushed along the polymer. In this case, 
the approximation of neglecting the switching on or off 
of the barrier (via electric field) holds well. As explained 
in Ref. for an electric field of 10 MV/cm (a typi- 

cal field for nanoscale field effect devices), the threshold 
gap frequency for PVDF corresponds to a temperature of 
roughly 5 K. Therefore, if the device operates at temper- 
atures below this range, these phonons will be effectively 
blocked from participating in heat transmission. Cou- 
pling PVDF to a grid of backgate electrodes that can be 
individually controled would effectively produce a mov- 
ing large barrier potential, as required by our pump. To 
evaluate the cooling capability of the pump, let us use 5 
K as an estimate for A (set by the threshold gap men- 
tioned above) and a velocity ratio vb/c = 1/10. Then, 
Sb ~ 0.4 K; it follows from Eq. ^ that if Th = 100 
mK, Tmin ~ 20 mK, while ff = 5 mK, Tnun « 40 /iK. 
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